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Abstract 

Consider a family of smooth potentials Vg, whicli, in the limit e ^ 0, become 
a singular hard-wall potential of a multi-dimensional billiard. We define auxil- 
iary billiard domains that asymptote, as e ^ to the original billiard, and provide 
asymptotic expansion of the smooth Hamiltonian solution in terms of these bil- 
liard approximations. The asymptotic expansion includes error estimates in the 

norm and an iteration scheme for improving this approximation. Applying this 
theory to smooth potentials which limit to the multi-dimensional close to ellip- 
soidal billiards, we predict when the separatrix splitting persists for various types 
of potentials. 

1 Introduction 

Imagine a point particle travelling freely (without friction) on a table, undergoing elas- 
tic collisions with the edges of the table. The table is just a bounded region of the plane. 
This model resembles a game of billiards, but it looks much simpler - we have only one 
ball, which is a dimensionless point particle. There is no friction and the table has no 
pockets. The shape of the table determines the nature of the motion (see f201 and 
references therein) - it can be ordered (integrable, e.g. in ellipsoidal tables), ergodic 
(e.g. in generic polygons), strongly mixing (in dispersing-Sinai tables or focusing- 
Bunimovich tables), or of a mixed nature for a general geometry with both concave 
and convex boundary components. A mechanical realization of this model in higher 
dimensions appears when one considers the motion of rigid (/-dimensional balls in a 
(/-dimensional box (c/ = 2 or 3): it corresponds to a billiard problem in a complicated 
n dimensional domain, where n = 2N x d (|29 jOJ ). 

Usually, in the physics context, this billiard description is used to model a more 
complicated flow by which a particle is moving approximately inertially, and then is re- 
flected by a steep potential. The reduction to the billiard problem simplifies the analysis 
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tremendously, often allowing to describe completely the dynamics in a given geometry. 
Numerous applications of this idea appear in the physics literature; It works as idealized 
model for the motion of charged particles in a steep potential, a model which is often 
used to examine the relation between classical and quantized systems (see 11711321 and 
references therein); This approximation was utilized to describe the dynamics of the 
motion of cold atoms in dark optical traps (see 1 18 1 and reference therein); This model 
has been suggested as a first step for substantiating the basic assumption of statistical 
mechanics - the ergodic hypothesis of Boltzmann (I21J. L29J . L30J . L31J .I33 1). The op- 
posite point of view may be taken when one is interested in studying numerically the 
hard wall system in a complicated geometry (e.g. apply ideas of [24'] to |25|) - then 
designing the "correct" limiting smooth Hamiltonian may simplify the complexity of 
the programming. 

For two-dimensional finite-range axis-symmetric potentials 1291 1221 l23l [11 1191 1151 
[0113^ it was shown that a modified billiard map may be defined, and several works 
have utilized this modified map to prove ergodicity of some configurations I29ll22ll23l 
[15','2 1, or to prove that other configurations may possess stability islands p |ll3l . The 
general problem of studying the limiting process of making a steep two-dimensional 
potential steeper up to the hard-wall limit can be approached in a variety of ways. In 
|24 1 approach based on generalized functions was proposed. In 1 35 1 we developed a 
different paradigm for studying this problem. We first formulated a set of conditions 
on general smooth steep potentials in two-dimensional domains (C smooth, not neces- 
sarily of finite range, nor axis-symmetric) which are sufficient for proving that regular 
reflections of the billiard flow and of the smooth flow are close in the topology. 
This statement, which may appear first as a mathematical exercise, is quite powerful. 
It allows to prove immediately the persistence of various kinds of billiard orbits in the 
smooth flows (see L35J and Theorem [S] in Section and to investigate the behav- 
ior near singular orbits (e.g. orbits which are tangent to the boundary) by combining 
several Poincare maps, see for example |27 36 9|. The first part of this paper (see 
Theorems lll2t is a generalization of this result to the multi-dimensional case. 

Thus, it appears that the Physicists approach, of approximating the smooth flow 
by a billiard has some mathematical justification. How good is this approximation? 
Can this approach be used to obtain an asymptotic expansion to the smooth solutions? 
The second part of this paper answers these questions. We propose an approximation 
scheme, with a constructive twist - we show that the best zero-order approximation 
should be a billiard map in a slightly distorted domain. We provide the scaling of 
the width of the corresponding boundary layer with the steepness parameter and with 
the number of derivatives one insists on approximating. Furthermore, the next order 
correction is explicitly found, supplying a modified billiard map (reminiscent of the 
shifted billiard map of 1291 1131 ') which may be further studied. We beUeve this part 
is the most significant part of the paper as it supplies a constructive tool to study the 
difference between the smooth flow and the billiard flow. 

Indeed, in the last part of this paper we demonstrate how these tools may be used 
to instantly extend novel results which were obtained for billiards to the steep po- 
tential setting; It is well known that the billiard map is integrable inside an ellipsoid 
1201 . Moreover, Birkhoff-Poritski conjecture claims that in 2 dimensions among all 
the convex smooth concave billiard tables only ellipses are integrable 1341 . In 1371 this 
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conjecture was generalized to higher dimensions. Delshams et al (" 111 11 . 1121 see refer- 
ences therein) studied the affect of small entire symmetric perturbations to the ellipsoid 
shape on the integrability. They proved that in some cases the separatrices of a simple 
periodic orbit split; Thus, they proved a local version of Birkhoff conjecture in the 2 
dimensional setting, and provided several non-integrable models in the n dimensional 
case. Here, we show that a simple combination of their results with ours, extends their 
result to the smooth case - namely it shows that the Hamiltonian flow, in a sufficiently 
steep potential which asymptotically vanishes in a shape which is a small perturba- 
tion of an ellipsoid, is chaotic. Furthermore, we quantify, for a given perturbation of 
the ellipsoidal shape, what "sufficiently steep" means for exponential, Gaussian and 
power-law potentials. 

These results may give the impression that the smooth flow and the billiard flow 
are indeed very similar, and so a Scientist's dream of greatly simplifying a complicated 
system is realized here. In the discussion we go back to this point - as usual dreams 
never materialize in full. 

The paper is ordered as follows; In Section |2] we define and describe the billiard 
flow and billiard map. In Section |3l we study the smooth Hamiltonian flow; we first 
prove that if the potential satisfies some natural conditions the smooth regular reflec- 
tions will limit smoothly to the billiard's regular reflections (Theorems lll2l i. Then, we 
define a natural Poincare section on which a generalized billiard map may be defined 
for the smooth flow. Next, we derive the correction term to the zeroth order billiard 
approximation (Theorem |3} and calculate it for three model potentials (exponential, 
Gaussian and power-law). We end this section by stating its immediate implication - 
a persistence theorem for various types of trajectories (Theorem|5|i. In Section0]we 
apply these results to the perturbed ellipsoidal billiard. We end the paper with a short 
summary and discussion. The appendices contain most of the proofs, whereas in the 
body of the paper we usually only indicate their main steps. 

2 Billiards in d dimensions 
2.1 The Billiard Flow 

Consider a billiard flow as the motion of a point mass in a compact domain D e or 
T''. Assume that the boundary dD consists of a finite number of C'+' smooth (r > 1) 
{d — 1) -dimensional submanifolds: 

aD = riur2u...ur„, i^\...n. (i) 

The boundaries of these submanifolds, when exist, form the corner set of dD\ 

r* = ariuar2u...uar„, i = \...n. (2) 

The moving particle has a position q £ D and a momentum vector p e M'' which are 
functions of time. If ^ e int{D), then the particle moves freely with the constant veloc- 
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ity according to the rule' : 

' q = p 



Equation (|3} is Hamiltonian with the Hamiltonian function (hereafter — {p,p)) 

H{q,p)^'-. (4) 

The particle moves at a constant speed and bounces of dD according to the usual elastic 
reflection law : the angle of incidence is equal to the angle of reflection. This means 
that the outgoing vector pout is related to the incoming vector /?,„ by 

Pout = Pin - 2{pi„,n{q))n{q), (5) 

where n{q) is the inward unit normal vector to the boundary dD at the point q, see IIIOI . 
To use the reflection rule (|5}, we need the normal vector n{q) to be defined, hence the 
rule cannot be applied at points ^ e F*, where such a vector fails to exist^. 

Definition 1. The domain D is called the configuration space of the billiard system. 

The phase space of the system isV ~Dx S''^^, where S''^ ' is a (li — 1 )-dimensional 
unit sphere (we set H = j) of velocity vectors. So the elements of fP are 

p = iq,p)- 

Denote the time t map of the bilhard flow as 

b,:po~ypt. (6) 

We do not consider reflections at the points of the corner set, so p, = btPo implies here 
that the distance between any point on the trajectory connecting with q, and the set 
r* is bounded away from zero. A point p G £P is called an inner point if q ^ dD and 
a collision point if q E dD \ F*. Obviously, if po and p, = b,po are inner points, then 
Pf depends continuously on po and t. If P( is a (non-tangent) collision point then the 
velocity vector undergoes a jump. Thus, in this case both bt-o and bt+o are defined. 
The map Ro = bt+ob^^Q is the reflection law (|5} (augmented by qour = qin). 

If the piece of trajectory that connects qo with qt does not have tangencies with the 
boundary, then p, depends C'-smoothly on po. It is well-known ( L30il . lt35J ) that the 
map bf loses smoothness at any point qo whose trajectory is tangent to the boundary 
at least once on the interval (0,f). Clearly a tangency may occur only if the boundary 
is concave in the direction of motion at the point of tangency. Consider hereafter only 
non-degenerate tangencies, namely assume that the curvature in the direction of motion 
does not vanish. 



' We assume that the particle has mass one (otherwise rescale time). 

^To be precise, one may define n{q) by continuity at points of F*, but this might give more than one 
normal vector n{q), hence the dynamics would be multiply defined for a genetic corner. We adopt a standard 
convention that the reflection is not defined at any £ F* . 
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Figure 1: Singularity near a tangent trajectory. For better visualization we present a 
slanted hyperplane which is divided into 2 parts: has a square-root singularity on the 
boundary between and Aj. 

Choose local coordinates q = {x,y) in such a way that the origin corresponds to 
the collision point, the j-axis is normal to the boundary and looking inside the bilUard 
region D, and the x-coordinates (x G M''^') correspond to the directions tangent to the 
boundary. If Q{x,y) = is the equation of the boundary in these coordinates, then 
2y(0,0) 7^ and QxiO,0) = 0. We choose the convention that 2y(0,0) > 0. Obviously, 
the tangent trajectory is characterized by the condition py = 0, where {px,Py) are the 
components of the momentum p. The vector Px = x indicates the direction of motion 
of the tangent trajectory. It is easy to check that the tangency is non-degenerate if and 
only if 

pIQxx{0,0)px>0. (7) 

Notice that if the billiard's boundary has saddle points (or if the bilUard is semi- 
dispersing), then there always exist directions for which this non-degeneracy assump- 
tion fails. On the other hand, if the boundary is strictly concave, then all tangencies are 
non-degenerate. 

Let x= (xi , . . . ,Xrf_i) with xi corresponding to the direction of motion (i.e. Px = 
(1 , 0, . . . , 0)). Then, the boundary surface near the point of non-degenerate tangency is 
described by the following equation: 

y = -ax] + 0{z^,xiz}, a > 0, 

where we denote z= (x2, . . . ,x,/_i). It is easy to see now that for a non-degenerate 
tangency, for a small z the map bx of the Une po = {xq = (— x/2, 0, . . . ,0), yo <0, pox = 
(1,0,..., 0), poy = 0) is given by 

Px = ((x/2,0, . . . ,0) + O(>'o), 2x^^+0(>'o), (1,0, . . . ,0) + O(>'o), 4,/^^+0{yo)). 
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As we see, the billiard flow looses smoothness indeed (it has a square-root singularity 
in the limit y —0) near the tangent trajectory. See Figure[2 

2.2 The Billiard map 

It is standard in dynamical system theory to reduce the study of flows to maps by 
constructing a cross-section. The latter is a hypersurface transverse to the flow. For 
the flow bt, such a hypersurface in phase space T can be naturally constructed with the 
help of the boundary of D, i.e. the natural cross-section S corresponds exactly to the 
collision points of the flow with the domain's boundary: 

5 = {p = {q,p) ev-.qedD, {p,n{q)) > 0}. (8) 

This is a {2d — 2) -dimensional submanifold in £P . Any trajectory of the flow bt crosses 
S every time it reflects at dD. This defines the Poincare map 

B-.S^S such that Bp = /Jt:o(p)+oP, (9) 

where 

To(p) = min{f > : /7,+oP S S}. 

Definition 2. The map B is called the billiard map. 

It is convenient to represent the billiard map as a composition of a free-flight and a 
reflection: 

B^RooFo, 

where the free-flight map is given by 

Fo{q,p) ^xo(p)-o(?:/'): (10) 
and the reflection law is given by 

Ro{q,p) = iq,P-2{p,n{q))n{q)). 

The billiard map B is a C — diffeomorphism at all points p e 5'\ E such that Bp € 5 \ E, 
where E is the singular set 

E = E,„„ge„„>s|jE£.„,.„e„ = eP: {p,n{q)) = 0}U {{q,p) <E P : q eT*}, (11) 

and B is C*^ at the non-degenerate tangent trajectories. 

3 Smooth Hamiltonian approximation 
3.1 Setup and Conditions on Potential 

Consider the family of Hamiltonian systems associated with: 

H^^+V{q;e), (12) 
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where the C^^' -smooth potential V{q;E) tends to zero inside a region D as e ^ 0, 
and it tends to infinity (or to a constant larger than the fixed considered energy level, 
say H — outside. Formally, the billiard flow in D may be expressed as a limiting 
Hamiltonian system of the form: 

Hh^'— + Vh{q), (13) 



where 



Let us formulate conditions under which this simplified billiard motion approximates 
the smooth Hamiltonian flow. In the two-dimensional case these conditions were in- 
troduced in 1351 . 

Condition I. For any compact region K <zD the potential V {q;£,) diminishes along 
with all its derivatives as £ 0: 

lmi||y(^;e)|,e^||c.+i =0. (15) 

The growth of the potential to infinity across the boundary needs to be treated more 
carefully. We assume that V is evaluated along the level sets of some finite function near 
the boundary. In other words, suppose, that in a neighborhood D of D\r* there exists 
a pattern function Q{q\€) : D ^ M} which is C^^ with respect to q and it depends 
continuously on e (in the C^^' -topology) at e > (so it has, along with all derivatives, 
a proper limit as e ^ 0). See Figure|2] Assume that away from F*: 

Condition Ila. The billiard boundary is composed of level surfaces ofQ{q;0): 

Q{q;£^ 0)|^Gr, = Qi = constant. (16) 

For each neighborhood of the boundary component F,- (so Q{q;t) is close to Qi), 
let us define a barrier function Wi{Q;z) : M} M}, which does not depend explicitly 
on q, and assume that: 

Condition lib. There exists a small neighborhood Ni of the surface F,- in which: 

V{q-z)\,eN. = W,(e(<?;e) -e,;e), (17) 

and 

Condition lie. Vy does not vanish in a finite neighborhood of the boundary sur- 
faces, thus: 

Ve|,e;v,^0 (18) 

and 

^W,(e-e,;e)^0. (19) 

Now, the rapid growth of the potential across the boundary may be described in 
terms of the barrier functions alone. Note that by ( I18> . the pattern function Q is mono- 
tonic across F,, so either Q > Qi corresponds to the points near F, inside D and Q < Qi 
corresponds to the outside, or vice versa. 
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Figure 2: Level sets of a pattern function Q{q;E). A bold line is a trajectory of the 
Hamiltonian flow near the boundary; a solid is a billiard trajectory. 
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Condition III. There exists a constant (may be infinite)'E > such that as e — > +0 
the barrier function increases from zero to "L across the boundary F,-; 



By Q could be considered as a function of W and e near the boundary: Q = 
Qi + Q, i(W';£). Condition IV states that for small £ a finite change in W corresponds 
to a small change in Q: 

Condition IV. As e — > +0, for any fixed W\ and W2 such that <W\ < W2 < E, 
for each boundary component F/, the function Q, i{W; e) tends to zero uniformly on the 
interval [Wi , W2] along with all its + 1 ) derivatives. 

Figure |3] shows the geometric interpretation of the pattern function and a typical 
dependence of the barrier function on Q and e. 

Note that the use of the pattern and barrier functions essentially reduces the d- 
dimensional Hamiltonian dynamics to a l-dimensional one, which allows for a direct 
asymptotic integration of the smooth problem. 

3.2 C° and - closeness Theorems 

Tlieorem 1. Let the potential V{q;e) in 11121 satisfy Conditions I-IV stated above. Let 
hf be the Hamiltonian flow defined by A121 on an energy surface H ~ H* < E, and 
bt be the billiard flow in D. Let po and pr — brPo be two inner phase points^. As- 
sume that on the time interval [0, T] the billiard trajectory of po has aflnite number of 
collisions, and all of them are either regular reflections or non-degenerate tangencies. 
Then hfp — ^^rp, uniformly for all p close to po and all t close to T. 

Tlieorem 2. In the conditions of Theorem^ further assume that the billiard trajectory 
of Po has no tangencies to the boundary on the time interval [0,T]. Then hf -^J^t in 
the C -topology in a small neighborhood of pa, and for all t close to T. 

The proof of the theorems is presented in the appendix and it follows closely the 
proof in 1 35 1. Informally, the logic behind Conditions I-IV is as follows. 

Condition I, obviously, implies that the particle moves with almost constant ve- 
locity (along a straight line) in the interior of D until it reaches a thin layer near the 
boundary where V runs from zero to large values (a smaller e corresponds to a thinner 
boundary layer). Note that the boundary layer can not be fully penetrated by the par- 
ticle. Indeed, as in all mechanical Hamiltonians, the energy level defines the region of 
allowed motion: for a fixed energy level // = //*< E , all trajectories stay in the region 
V(^;e) <H* . It follows from Condition III that for any such H* , the region of allowed 
motion approaches D as £ ^ 0. Thus, by Condition III, if the particle enters the layer 
near a boundary surface (note that points from F* are not considered in this paper), it 
has, in principle, two possibilities. First, it may be reflected and then exits the bound- 
ary layer near the point it entered. The other possibility, which we want to avoid, is 
that the particle sticks to the boundary and travels along it far from the entrance point. 
Condition IV guarantees that if the reflection is regular, or in case of non-degenerate 

^Hereafter, T always denotes ?i finite number. 




Q inside D 
Q outside D 



(20) 
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tangency, the travel distance along the boundary vanishes asymptotically with e. The 
case of degenerate tangencies, which are unavoidable in the higher dimensional case 
if the boundary has directional curvatures of opposite signs (namely saddle points), is 
not studied here. Once we know that the time spent by the particle near the boundary 
is small, we can see that Condition II guarantees that the reflection will be of the right 
character, namely the smooth reflection is C^-close to that of the billiard. Indeed, Con- 
dition II implies that the reaction force is normal to the boundary, hence, as the time 
of collision is small and the position of the particle does not change much during this 
time, the direction of the force stays nearly constant during the collision. Thus, only the 
normal component of the momentum is changing sign while the tangent components 
are nearly preserved. Computations along these lines provide a proof of Theorem[2 

Proving Theorem|2] i.e. the C^-closeness, makes a substantial use of Condition IV. 
Let us explain in more detail the difference between the C'^ and topologies in this 
context. Take the same initial condition {qo,po) for a billiard orbit and for an orbit of 
the Hamiltonian system J12t (the Hamiltonian orbit will be called the smooth orbit). 
Consider a time interval t for which the billiard orbit collides with the boundary only 
once. In these notations (p,„ is the angle between po (the momentum at the point ^o) 
and the normal to the boundary at the collision point, (po„, is the angle between p, (the 
velocity vector at the point q,) and the normal. Define the incidence and reflection 
angles ((p,„(e) and (Po„,(e)) for the smooth trajectory in the same way. Theorem ^ 
implies the correct reflection law for smooth trajectories: 

(p,>,(e) + (p„„f(e) «0 (21) 

for sufficiently small 8. However, (p,„ + (p„j„ is a function of the initial conditions, so a 
non-trivial question is when it is close to zero along with all its derivatives. In Theorem 
121 we prove that Condition IV is sufficient for guaranteeing the correct reflection law 
in the C' -topology in the case of non-tangent collision (near tangent trajectories the 
derivatives of the smooth flow cannot converge to those of the billiard because the 
billiard flow is singular there, see Figure[l). 

Hereafter, we will fix the energy level of the Hamiltonian flow to H* = j. Notice 
that the analysis may be applied to systems with steep potentials which do not depend 
explicitly on e (or do not degenerate as e — > 0) in the limit of sufficiently high energy: 
the reduction to the setting ( I12t which we consider here may be achieved by a scaling 
of time. 

3.3 Asymptotic for a regular reflection 

It follows from the proof of Theorem |2] that the behavior of smooth trajectories close 
to billiard trajectories of regular reflections can be described by an analogue of the 
billiard map. More precisely, one can construct a cross-section in phase space of 
the Hamiltonian flow, close to the "natural" cross-section S where the billiard map B 
is defined; the trajectories of the Hamiltonian flow which are close to regular billiard 
trajectories define the Poincare map on S^, and this map is C^-close to B. Let us explain 
this in more details. 

It is convenient to consider an auxiliary billiard in the modified domain D^, defined 
as follows. For each boundary surface F,, take any V, (e) +0 such that the function 
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Figure 4: Free flight between boundaries F? and F^. A smooth trajectory is marked by 
a bold hne and an auxiliary billiard trajectory is marked by a solid line. 

(inverse barrier) Q, ,(W;e) tends to zero along with all its derivatives, uniformly for 
5 > > V,-. We will use the notation 

m\''\v-z)= sup \Q}'\w■,^)\. (22) 
V < W < i 
< / < r+ 1 

Condition IV implies that M approaches zero as e ^ for any fixed v > 0, hence 
the same holds true for any sufficiently slowly tending to zero v(e), i.e. the required 
V;(e) exist. Let 'n,(e) = Q, ,(v,;e) and consider the billiard in the domain which is 
bounded by the surfaces F? : Q{q;t) = +r|,(e). See Figure|4] Recall that the bound- 
aries F, of the original billiard table D are level sets Q{q;Q) — Qi, and that r|,(e) 
by construction, so the new billiard is close to the original one. In particular, for regu- 
lar reflections, the billiard map of the auxiliary billiard tends to the original billiard 
map B along with all its derivatives. It is established in the proof of Theorem |2] that 
for any choice of v,'s tending to zero, the condition q e dD^ defines a cross-section 
in the phase space of the smooth Hamiltonian flow; Trajectories which are close to 
the billiard trajectories of regular reflection, i.e. those which intersect dD^ at an angle 
bounded away from zero, define the map 

: {q £ dD^,p looking inwards D^) {q e dD^,p looking outwards D^), 

namely 

F\q,p)=h%{q,p) (23) 
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Figure 5: Reflection from the boundary F,-. A smooth trajectory is marked by a bold 
line. An auxiliary billiard trajectory only changes its direction according to the law 
(|23. 

and this map is close to the free-flight map (see Section l2!2l of the billiard in D^: 

F^{q,p)=b,._,iq,p) (24) 

where X^{q,p) is the time the smooth Hamiltonian orbit of {q,p) needs to reach dD^, 
and '^{q,p) denotes the same for the billiard orbit. Note that we cannot claim the 
closeness of the time T maps for the smooth Hamiltonian and billiard flows everywhere 
in D^, still we claim that the maps ( I23> and (I24> are close; we will return to this later. 

Outside D^, the overall effect of the motion of smooth orbits is close to that of 
a billiard reflection. Namely, as it is proved in Theorem |2] once v, is chosen such 

(r) 

that Mj (v,-,e) 0, the smooth trajectories which enter the region W,(2;e) > V,- at a 
bounded away from zero angle to the boundary, spend in this region a small interval 
of time (denoted by T^(<7,„, /?,„)) after which they return to the boundary W,(2;e) = V,- 
(namely to Q{q;E) = 2, +T|,(e)). Thus, these orbits define the map 

: [qin e dD^,pi„ looking outwards D^) — > e dD^ ,paut looking inwards D^). 

It follows from the proof of Theorem |2l that the map R'^ is close to the standard 
reflection law R% from the boundary dD'^ : 

Rl{q^p) = {q,P ~ 2n{q) {n{q),p)), (25) 

where n{q) is the unit normal vector to the boundary dD^ at the point q. See Figure 
|5] Note that the smooth reflection law R^ corresponds to a non-zero (though small) 
collision time x^{q,p), unlike the billiard reflection which happens instantaneously. 
Summarizing, from the proof of Theorem|2]we extract that on the cross-section 

s, = {p = {q,p) -.qedD^ {pAq)) > 0} (26) 
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the Poincare map 

^^=R^o (27) 

is defined for the smooth Hamiltonian flow (for regular orbits - orbits which intersect 
dD^ at an angle bounded away from zero), and this map is C^-close to the billiard map 

= RloF^. As the billiard map is close to the original billiard map B, we obtain the 
closeness of the Poincare map (S>^ to B as well. However, when developing asymptotic 
expansions for <I>^, it is convenient to use the map B'^ (rather than B) as the zeroth order 
approximation for <P^. Then, the next term in the asymptotic may be explicitly found 
(see below) and the whole asymptotic expansion may be similarly developed. 

We start with the estimates for the "free flight" segment of the motion, i.e. for 
the smooth Hamiltonian trajectories inside D^. For every boundary surface F,-, choose 
some 5,- (e) — > such that the surfaces Q{q;e) ^ Qi + 5, (e) bound the region inside 

in which the potential V tends to zero uniformly along with all its derivatives. See 
Figure |6l Let 

mW(5;e)= sup \\d'v{q;z)\\. (28) 

1 </<r+l 

According to Condition 1, m approaches zero as e ^ for any fixed 5 of the appropriate 
signs, therefore the same holds true for any choice of sufficiently slowly tending to 
zero 5, (e). As m''^' 0, it follows that within the flow of the smooth Hamiltonian 
trajectories is C^-close to the free flight, i.e. to the billiard flow. In other words, the 
time T map h^{q,p) = {qz,Pi) of the smooth flow in is 0^,^ (m'''')-close to the time 
1 map of the billiard flow 

b.{q.p)={^\P'). (29) 
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Note that on the boundary of we have, by construction, qJ- (IV; e) — > 0, i.e. W'j (2; e) ^ 
oo, while on the boundary of D^^^^ we have W/(2;e) 0. Thus, we have a boundary 
layer D^\Dl^f of a non-zero width |5,(e) — r|,(£)| in which the gradient of the poten- 
tial rapidly decreases. The speed with which the value of Q{q{t);e) changes within 
this boundary layer is bounded away from zero (see the proof of Theorem |2j, so the 
time the orbit needs to penetrate it is (9(5,). Within this boundary layer the time x map 
{q,p) I— > (q^jP-c) of the smooth flow is not necessarily close to the time T map of the bil- 
liard flow i29\ . However, it is shown in the proof of Theorem|3] that the maps from one 
surface Q = const to any other such surface within the boundary layer are C^-close for 
the two flows. This, obviously, implies the closeness of the maps and (because 
the corresponding cross-section is the surface of the kind Q = const indeed). 

In Appendix 17. 2l we show that by an appropriate change of coordinates in each of 
the three regions we consider (inside D^^,, in D'^\Dl^f, and outside D'^), the equations of 
motion may be written as differential equations integrated over a finite interval with a 
right hand side which tends to zero in the C^-topology as e ^ 0. Thus, not only do we 
obtain error estimates for the zeroth order approximation, we also find a method for ob- 
taining higher order corrections using Picard iterations; The asymptotic behavior of the 
right hand side of the equations leads to a contractivity constant which asymptotically 
vanishes and thus the Picard iteration scheme provides asymptotic for the solutions 
(each new iteration provides a better asymptotic). In this way we prove in appendix l7.3l 
the following 

Lemma 1. Let q be an inner point of D, and p be such that the first hit of the billiard 
orbit of {q,p) with the boundary F,- is non-tangent. Then, the orbit of the smooth flow 
hits the cross-section {q e F?} — {Q{q;z) — 2i the point (q'x,/'x) such that 

q + px + 0^.,{m'^''^ +Vi) 
q + pT + Vy {q + ps- e) {s -%)ds + 0^^_, ( (m('-) + V^^) , 

P + 0„.(mW+V,) ^^^^ 
P - Io^V{q + ps;e)ds + 0^^_, ((mW + v,)^), 

where 1 — l^{q,p) denotes the travel time to the boundary ofD^ (so 2(^x;£) = Qi + 
X\q,p) = T^(^,/?) + 0^,(mW+V,) 

= + (31) 

where VQ is evaluated at the (auxiliary) billiard collision point q + pl%{p,q), and 
'l%{p,q) is the time the billiard orbit of [q, p) needs to reach F^. 

Now, let us estimate the free-flight map F^ of the Hamiltonian flow. If ^ e F^ and 
{p,n{q)) is positive and bounded away from zero, and if the straight line issued from q 
in the direction of p first intersects dD^ (say, the surface F?) transversely as well (in our 
notations this can be expressed as the condition that {p,n{q + pll{q,p))) is negative 
and bounded away from zero), then the orbits of the Hamiltonian flow define the map 
F^ from a small neighborhood of {q,p) on the cross-section {q e F^} in phase space 



?x = 



14 



into a small neighborhood of the point {q + pzl{q,p),p) on the cross-section {q G F?}. 
See Figure 0] Take an inner point {qi,pi) on the smooth Hamiltonian trajectory of 
{q,p). By construction (see dl^). 

1^iq,p)=Z^{qu-pi)+Z^iqi,pi), 

{q,p)=h%,^^^_p^^{qi,pi) 

and 

As qi is bounded away from the billiard boundary, we can plug ( 13 0> and ( I31> in these 
relations, which gives us the following 

Lemma 2. Near the point {q,p) under consideration, the free flight map : {q,p) ^ 
{q^,p^) for the smooth Hamiltonian flow is O^r {m'''^ + Vi + V j)-close to the free flight 
map of the billiard in and is given by 

q^ =q + pX^ + VV{q + ps■^){s ~ X^)ds + O +Vi + Vjf), 
P,e=P- /(f '7V{q + ps;E)ds + O ((m('-) + V, + Vy)^). 

The flight time X^{q,p) is O^y (m^''^ + V; +V j)-close l^{p,q) and is uniquely deflned by 
the condition Q^q^^'.^) — Qi+y\i{'£) (cf.UTl): 

9, X 9, . {VQj^°VV(q + ps;e)(zl~s)ds) 
X'iq,p) = Xl{q,p) + ^-^^ "^^^J^^ " LJ. + o^,_, ((m(') +V, +V,)'), 

(33) 

where VQ is taken at the billiard collision point q + pl%{p,q) andi^^{p^q) corresponds 
to the free flight travel time: Q{q + pll{p,q)',E) ~ Qi + r|,(e). 

This could be written as 

F' = F! + O^., (mW + V, + V,-) = F^+Fl + 0^^_, ((m*'") + V,- + V;)'), 

where Ff = O^., {m^''^ + V; + V;) and F^ is defined bvl24l 

Note that the above estimates hold true for any choice of 5,'s such that m^''^ 0. 
Therefore, one may take 5,'s tending to zero as slow as needed in order to ensure as 
good estimates as possible for the error terms in ( I32> . (I33> . 

Next we estimate the reflection law for the smooth orbit. Consider a point <? S F? 
and let the momentum p be directed outside D^, at a bounded from zero angle with Ff. 
As we explained, the smooth trajectory of {q,p) spends a small time T^{q,p) outside 

and then returns to Ff with the momentum directed strictly inside D^. Let py and 
Px denote the components of momentum, respectively, normal and tangential to the 
boundary F^ at the point q: 

Py = {n{q),p), Px= P- Pyn{q). (34) 
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We assume that the unit normal n{q) is oriented inside D^, so py < at the initial point. 
Denote by Qy{q;£) the derivative of Q in the direction of n{q): 

Qyiq;e):={VQ{q;e),niq)), 

let K{q;e) denote the derivative of n{q) in the directions tangent to Ff, and let l{q;e) 
denote the derivative of n{q) in the direction of n{q). Obviously, Qy is a scalar, K is 
a matrix and / is a vector tangent to Ff at the point q. Note that Qy ^ by virtue of 
Condition lie. Define the integrals: 



I,^I,{q,p)^2j,'^(l',C-J^;e)ds 



(35) 



and the vector J: 

J, X h{q-,p) 
J{q,p) - 



l{q;e)+Ii{q,p)K{q;e)p, 



/Qy{q;e). (36) 



Py 

Notice that 7 is a vector tangent to Ff at the point q and that by i22i . 

h ,2 = O^, (Mf ) ) , y ^ o^^_ , (m|'"' ) . (37) 
In Appendix l7.3l we prove the following 
Lemma 3. For the smooth Hamiltonian flow, the collision time is estimated as 

zUq,p) = = -^-L-^I,iq,p) + 0^^_, mh\ (38) 

The reflection map : {q,p) i— s- {q,p) is given by: 
q = q + 0^.,{M^h =q + p,zl{q,p)+0^,._, {{M^^ f), 

p = p- 2n{q)p, + O^r (M,*'*) =p~ 2n{q)p, - PyJiq,p) - n{q) {p_„J{q,p)) + 0^^_, ((M?'' 

(39) 

As we see from this lemma (see also i3H ). 

= Rl + O,, (m('-)) = Rl+R\ + 0^^_, ((Mf))2), 

where R^ = O^^ ^ {M\ ) and R\ is defined by|25l Thus, the smooth reflection law is 

O^r (M|'^)-close to the billiard reflection law ( l25t . 
Combining the above lemmas we establish: 

Theorem 3. Let the potential V{q;€) satisfy Conditions I-IV, and choose 5,- 's and v, 's 
such that 5,(£),V,(e),m(''^(8),M?'^'' (e) — > as £ — s- 0. Then, on the cross-section (see 
i26l ) near orbits of a regular reflection'^, for all sufficiently small e the Poincare map 



■*that is, given any constant O 0, near tlie points {q,p) G Se such tliat {n{q),p) > C and | {n{q),p) \ > C 
where {q,p) =B^(q,p) 
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Figure 7; Billiard map B (red). Billiard map of the auxiliary billiard (green). 
Poincare map for the smooth Hamiltonian flow (blue). The first approximation 
of 4>'=B'= + 4>f (violet). 

of the smooth Hamiltonian flow is defined, and it is 0{m^^^ +V -\- M^''^)-close in 
the C -topology to the billiard map B^ = R^o in the auxiliary billiard table ( see 
Figure^. Furthermore, 

^^^R^oF^^B^ + 0^,{m^'''^+V+M^'^) = {Rl+R\)o{F!+Fl) + 0^,._,{{m^''^+V+M^''^f) 

=: +4>f + 0^^_^ ((m*'") +V+M^''^f) (40) 

(where V = max,V/, M^''^ = max,M|'^',4>j = ^c'-' +V +M^''^), and the first order 
corrections F^ and R\ are explicitly calculated in Lemmas^and]3\. 

Theorem 4. Given a finite T and a regular billiard trajectory in [0, T], the time t 
map of the smooth Hamiltonian fiow and of the corresponding auxiliary billiard are 
0(y + m^''^ + M^'''^)-close in the -topology for all t G T\Tr, where Tr is the finite 
collection of impact intervals each of them of length 0(|5| +M'''^). 

3.3.1 Error estimates for some model potentials 

Now we can estimate the deviation of the smooth Hamiltonian trajectories from the reg- 
ular (non-tangent, non-corner) billiard ones for various concrete potentials V{q;e). To 
make a general estimate possible, we have to assume that the behavior of the potential 
near the boundary dominates the estimate; We say that V{q;e) is boundary dominated, 
if y(g';8) and its derivatives are smaller in the interior of D?,^ (i.e. in the region bounded 
by the surfaces Q{q;e) = Qi + 5,(£)) than on the boundary of this domain. This means 



17 



that for boundary dominated potentials 

mW(5;e)= sup ||a'y(^;£)|| - sup ||a'y(^;e)||. (41) 

1 </<'-+! l</<r+l 

By the definition of the pattern function Q, near a given boundary F, 



V[q-Z] 



= W^(e(^;e)-e,;8) 



e=e,+8, 



Since Q{q;t) is bounded with its derivatives, we conclude that there exists a constant 
C such that 

m(')(5;e) =Cmax max |w/'V5,;£)|. (42) 

;■ l</<r+l 

Thus, for boundary dominated potentials, one can estimate the differences /jf — b, and 
cj)e _ in terms of the barrier functions alone. 

The corresponding estimates given by Theorems|3and|3]hold true for every choice 
of V and 5 such that 5(e),v(e),mr(5(e);£),M,-(v(£);£) ^ as £ ^ (for simplicity of 
notation we assume hereafter that the barrier function W is the same for all boundary 
surfaces F,-, and thus suppress the dependence on /). To obtain the best estimates, we 
have to find v(£) and 8(£) which minimize the expression v + M'^'"' (v; £) + w'^'"' (5;e). 
In this way, we first find v(£) which minimizes v +m('"' (v;£)- As m'*"' is a decreasing 
function of v (see \22V ). the sought v(£) solves the equation 

v = mW(v;£). (43) 

After V is determined, we may try to make 5(8) go to zero so slow that the correspond- 
ing value of m'^'"' (see (I42» will be asymptotically equal to v(£). Once succeeded, we 
may conclude that v(£) given by (I43> estimates the deviation between regular billiard 
and smooth trajectories. Notice that the significance of v(£) is three-folded; First, it 
determines the optimal auxiliary billiard which supplies the best approximation to the 
smooth Hamiltonian flow (see Lemma|3. Second, it estimates the accuracy of this ap- 
proximation. Third, it determines, via the relation m'''(5) ~ V, the width |5(£)| +v(£) 
of the boundary layer in which the billiard and the Hamiltonian flows are not close 
(Theoreml^. Let us proceed to examples. 

Proposition 1. Consider the boundary dominated potential V{q\z) corresponding to 
the barrier function W{Q) — e^~ for small Q. Then, near regular billiard trajectories, 
the smooth Hamiltonian flow is 0{ ''^\/s,)-close in the -topology to the billiard flow 
within the auxiliary billiard defined by the level set Q{q;e) = ri(£) = C>(£ln£). The 
corresponding Poincare map <i>^ is O^^r ( '^y^)-close to the auxiliary billiard map B^. 
The impact intervals lengths are 0{ '^-^/e). 

Proof Since W(')(g;£) = (-£)"'e"e, we obtain fliatmW(5;£) = 0(£"(''+''e"t) (since 
the potential is boundary dominated, we may use ( I42» . The inverse to W{Q;e) is 



18 



given by (1{W-e) = -elnW, so Q,(')(H^;e) = (-l)'(Z - l)!eW^-', and MW(v,e) = 
0(ev"(''+')) (see Plugging this in (|43}, we find 

v(e) = (44) 

By choosing 5(e) = — (r + 1 + ^)elne, we obtain m(''^(5,e) ^ v(e), so for v given by 
J44t we have that v+M'^'') = (9(v), and the proposition now follows immediately 

from Theorems|3]and|3(the value of ri(£) — (9(elne) is given by r| = Q,(v;e)). □ 

Proposition 2. Let the boundary dominated potential V{q;z) correspond to the bar- 

rier function W {Q) = e e for small Q. Then, near the regular billiard trajectories, 
the smooth Hamiltonian flow is (9(v(e)) — 0{ -close in the -topology to 

the billiard flow within the auxiliary billiard defined by the level set Q{q;e) = r|(e) = 
(9(-\/e|lne|). The corresponding Poincare map is O^, {v{E))-close to the auxiliary 
billiard map B^. The impact intervals are of the length 0{y{€)). 

Proof Itiseasy tosee thatVK(')(g;e) =(9((f )'e"^) forg> hence mW(5;e) = 
0{{\y+^e-^). From Q_{W\z) = V-elnlVwe obtain M('')(v;e) = C'(y^v-(''+i)). 
Plugging this in j43> . we indeed find 



as required. By choosing 5(e) ^ y— 5(r+ 1 + ■^)elne, we obtain m('''(5;£) ^v(e), 
so the rest follows directly from Theorems|3]and|3 □ 

Proposition 3. Let the boundary dominated potential V{q;€) correspond to the bar- 
rier function W{Q) — (^)"- Then, near the regular billiard trajectories, the smooth 

Hamiltonian flow is 0(v(e)) = 0{ '^^"^k/^)-close in the -topology to the billiard flow 
within the auxiliary billiard defined by the level set Q{q;e) = ri(£) — 0(v'^^^). The 
corresponding Poincare map is Oi^, {v{€))-close to the auxiliary billiard map B^. 

a(r+2) 

The impact intervals are (9(v(e)) when a > 1, and 0{v{e) o-+r+i ) when a < 1. 

Proof As above, using w(')(g;e) = 0{-^) we obtain that m('')(5;e) = 0(^1^), 
and since Q,(VK;e) = -j^, we find Q,(''(W;£) = 0{ ^,li/a ) and thus mW(v;£) = 
0{ ,-+i+i/n ). It follows that v(e) = 0{ solves v = MW(v;e). Now r|(e) = 

a(r+2) . . 

= ^ = 0(v''+^). By taking 5(e) = vnTFa, we ensure that m('"'(5,e) ~ v(e). 
The length of impact intervals is now given by 0(v + 5). □ 

Note that the asymptotic for the deviation of the smooth trajectories from the bil- 
liard ones and for the length of the impact intervals depend strongly on r, i.e. on the 
number of derivatives (with respect to initial conditions) which we want to control. 
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Figure 8: P'' is a billiard periodic orbit (solid). is a periodic orbit of the smooth 
Hamiltonian flow (bold). 

3.4 Persistence of periodic and homoclinic orbits 

The closeness of the billiard and smooth flows after one reflection leads, using standard 
results, to persistence of regular periodic and homoclinic orbits. For completeness we 
state these results explicitly: 

Theorem 5. Consider a Hamiltonian system with a potential V{q,€) satisfying Con- 
dition I-IV in a billiard table D. Let P^{t) denote a non-parabolic, non-singular pe- 
riodic orbit of a period T for the billiard flow. Then, for any choice o/v(e),5(£) 
such that v(e),5(8),m'-'^(e),M(^^(£) oi £ — > 0, the smooth Hamiltonian flow has 
a uniquely defined periodic orbit P^{t) of period — T -\- Oiy + m^^' +M^^''), which 
stays 0{v + m^^^ + M^^'^)-close to P^ for all t outside of collision intervals (finitely 
many of them in a period) of length 0(|5| +M^^-'). Away from the collision intervals, 
the local Poincare map nearP^ is O^, (v + m^''^ +M^''^)-close to the local Poincare map 
near P^. In particular, if P^ is hyperbolic, then P^ is also hyperbolic and, inside D^, 
the stable and unstable manifolds of P^ approximate O^r (v + m^''^ + M^''^)-closely the 
stable and unstable manifolds ofP^ on any compact, forward-invariant or, respectively, 
backward-invariant piece bounded away from the singularity set in the billiard' s phase 
space; furthermore, any transverse regular homoclinic orbit to P^ is, for sufficiently 
small e, inherited by P^ as well. 

As P'^ is a regular periodic orbit, i.e. it makes only regular reflections from the 
boundary (a finite number of them on the period), it follows from Theorem |3] that a 
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Poincare map for the smooth HamOtonian flow near is Oiy + m^'' +M('')-close in 
topology to the Poincare map of the auxiliary billiard D^, while the latter is (9(r|(e))- 
close to the Poincare map for the original billiard D. Moreover, from \22l it follows that 
ri(e) <mW <m(i) and we can conclude that a Poincare map for the smooth Hamilto- 
nian flow near is (9(v + m^'^ +M('')-close in C' topology to the Poincare map for 
the original billiard D. Since, by assumption, P''{t) is non-parabolic, the corresponding 
fixed point of the Poincare map persists for sufficiently small £ in virtue of the implicit 
function theorem (the closeness of the corresponding continuous-time orbits is given 
by Theorem^}. The continuous dependence of the invariant manifolds of e in the hy- 
perbolic case follows from the continuous dependence of the Poincare map on e at 
all e > (TheoremO, and implies the persistence of transverse homoclinics immedi- 
ately. Indeed, the formulation regarding the closeness of compact pieces of the global 
stable and unstable manifolds may be easily verified by applying finite time extensions 
of the local stable and unstable manifolds. Note that similar persistence result holds 
true for topologically transverse homoclinic orbits. 

More generally, one may claim (by the shadowing lemma) the persistence of com- 
pact uniformly hyperbolic sets composed of regular billiard orbits. Note that the ac- 
curacy of the approximation of smooth orbits (periodic and aperiodic) by the billiard 
ones, does not depend on the orbit (e.g. is independent of its period) and is given by 
the maximal deviation for each reflection (times a constant). This holds true for any 
compact set of regular orbits of a strictly dispersing billiard flow (since such billiards 
are uniformly hyperbolic); see for example a nice application by Chen 1 9 1 . 

In some cases, to establish the existence of transverse or topologically transverse 
homoclinic orbits in a family of billiard flows bt (y) in Dy, one uses higher dimen- 
sional generalizations of the Poincare-Melnikov integral (see Section|3i. In particular, 
with the near integrable setting, the "splitting distance" between the manifolds near 
the transverse homoclinic orbit may be proportional to an unfolding parameter y. The 
above theorem impHes that if eo = eo(y) is chosen sothat v(£o,y)+»j'''(5(£o,y);eo,y) + 
M'''(v(£o,y);eo,y)) = o(y) and eo(y) ^ as y^ then, for sufficiently smaU y, trans- 
verse homoclinic orbits appear in the smooth flow for all e S (0,eo(y)). In the next 
section we use this remark and 1 11 1 to establish that transverse homoclinic orbits ap- 
pear in families of smooth billiard potentials which limit to the ellipsoidal billiard. 

4 Application to ellipsoidal billiards with potential 

Consider the billiard motion in an ellipsoid 

D = {qeW:{q,A-^q)<l}, (45) 

A — diag(t/i , . . . , dn) di > . . . > d„>Q. 

The ellipsoid is called generic if all the above inequalities are strict. A well known 
result of Birkhoff 1 3 1 is that the billiard motion in an ellipsoid is integrable, and the 
mathematical theory which may be invoked to describe and generalize this result is 
still under development - see Radnovic 1 16| and references therein. Delshams et al 
fl II and recently Bolotin et al |4| (see also references therein) investigate when small 
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non-quadratic symmetric perturbations to the ellipsoidal shape change the integrability 
property. In this series of works the authors prove the persistence of some symmetric 
homoclinic orbits, and for specific cases they prove that these orbits are transverse 
homoclinic orbits of the perturbed billiard, thus proving that integrability is destroyed. 
Here, we show that using the machinery we developed we can immediately extend 
their work to the smooth billiard-potential case (notice that in some results are 
extended to billiards with a C^-small Hamiltonian perturbation in the domain's interior, 
however the billiard potentials which we consider do not fall into this category - near 
the boundary they correspond to a large perturbation even in the C'-norm). We will 
first explain the relevant main results of Delshams et al, then supply the corresponding 
proposition for the smooth case (consequences of Theorem|2] or more specifically of 
Theorem O and then the corresponding quantitative estimates for specific potentials 
(which follows from Propositions lll3> . 

4.1 The billiard in a perturbed ellipsoid 

Consider the simplest unstable periodic orbit in an ellipsoidal billiard - the orbit along 
the diameter of the ellipsoid joining the vertices {—di,Q, ... ,0) and (iii,0, . . . ,0). De- 
note the set formed by the two-periodic points associated with the diameter by 

^' = {P+,P-} P± = {q±,P^} q±^i±diA---,0) p± = (±l,0,...,0). (46) 

These points correspond to isolated two-periodic hyperbolic orbits of the Billiard map 
B and the corresponding periodic orbit P/' = ^r(p+) of the billiard flow. The n — 1- 
dimensional («-dimensional for the flow) stable and unstable manifolds of this periodic 
orbit coincide; In 2-dimensions there are 4 separatrices connecting {p+,p_} whereas 
the topology of the separatrices in the higher dimensional case is non-trivial - it is well 
described by CW complexes for the 3 dimensional case and by hierarchal structure of 
separatrix submanifolds in the higher dimensional case (see im i. 

Of specific interest are the symmetric homoclinic orbits - it is established in II II 
that in the generic 2 dimensional case there are exactly 4 homoclinic orbits which are 
.X— symmetric (symmetric, in the configuration space, to reflections about the jc-axis) 
and 4 which are y— symmetric. In the generic 3 dimensional case, in addition to the 
16 planar symmetric orbits (8 in each of the symmetry planes- xy and xz) there are 16 
additional symmetric spatial orbits - 8 are symmetric with respect to reflection about 
the xz plane and 8 are y axial. In the n dimensional case there are 2"+' spatial symmetric 
orbits. ^ 

Denote by p^^hom _ |p6-hom| these symmetric homoclinic orbits 

of the billiard map in the ellipsoid, so p^-ho™ = spb-hom pb-hom ^ ^^^pb-hom-^ 
denotes the corresponding continuous orbit of the billiard flow. Given a q such that 
< i; <C c/n, define the local cross-sections of the billiard map by: 

{{q,p)\q^'dD,qi +di <q, l-pi< q}, 
= {iq,p)\qedD,di-qi < q, pi + l<q}, 

so, in particular, p± £ and C S, where S is the natural cross-section on which 
the billiard map is defined (see Section IT^ . It follows that only a finite number of 
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Figure 9: Billiard trajectory giving rise to a y-symmetric homoclinic orbit (blue). 



points in do not fall into E*, and that for any given geometry there exist a finite 

q such that p^'-hom^jpi-homp j;±j. _^ fQj. ^jj jjjg syimnetric orbits. See Figure |9l 
Thus, it is possible to choose p*'-hom ^ j^^^j cross-section E" such that f>*-hom ^ 
E** C {^yiE^ U E^}}. Notice that for the ellipsoid all the reflections are regular, and 
furthermore, for the symmetric homoclinic orbits, if d\ is finite and d„ is positive then 
all the reflection angles of p^^ho™ are strictly bounded away from %/2. 
Now, consider a symmetric perturbation of the ellipsoid Q of the form: 

D,^{qeW: {q,A-'q} < 1 +YS(|, • • • , |)}, (47) 

where the hypersurface Dy £ M"+' is symmetric with regard to all the coordinate 
axis of W and the function E : M" ^ M is either a general entire function, such that 
E(0, ... ,0) = or of a specific form (e.g. quadratic). By using symmetry arguments, 
Delshams et al prove that for generic billiard the above mentioned symmetric ho- 
moclinic orbits persist under such symmetric perturbations. Furthermore, analyzing 
the asymptotic properties of the symplectic discrete version of the Poincare-Melnikov 
potential (the high dimensional analog of the integral), they prove that for sufficiently 
small perturbations (small y) the n-dimensional symmetric homoclinic orbits are trans- 
verse in the following four cases: 

1. In two-dimensions, for narrow ellipses (Pi = -4 ^ 1), for any analytic small 
enough symmetric perturbation. 

J> 2 4 

2. In two-dimensions, in the non-circular case (pi ^ 1), for ^) = 
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3. In the three-dimensional case, for nearly flat elUpses (|32 = "i- ^ 1)^ for perturba- 

9 2 2 2 2 2 

tions of the form: ^j, ^i) = ^Ri^ijri) where is a generic polynomial 

1 z 3 3 2 3 

(or of some specific list). 

^2 

4. In the three-dimensional case, for nearly oblate elhpses (pi = -4 ~ 1), for the 
perhirbationS(|,J,|) = ||. 

To establish these results, the Poincare-Melnikov potential is calculated for each 
of these cases, and it is shown that it has non-degenerate critical points at the corre- 
sponding symmetric trajectories. It follows that p^-hom-y pgj-gjgts and the change in 
the splitting distance between the separatrices 'W " and W near Pq ^ is propor- 
tional to y, the perturbation amplitude, so that near Pq ^ at the local cross-section 

d {W^ ,w;')^M{r)y+0{f). (48) 

where r E R"^^ denotes some parametrization along W and M{r) (the gradient of the 
Poincare-Melnikov potential) has simple zeroes at the parameter values corresponding 
to any of the spatial symmetric homoclinic orbits Pq^^""^. 



4.2 Smooth Potential in a near ellipsoidal region 

Let us now consider a two parameter family of smooth potentials V{q;'Y, e) which limit, 
as e ^ to the billiard flow in the perturbed ellipsoid family Dy, namely, consider the 
family of Hamiltonian flows: 

H{e,y) = '— + V{q;y,e). (49) 

where y(^;y,e) satisfies conditions I-IV for all y values. In the four cases mentioned 
above, the flow limits, as e ^ 0, to an integrable billiard motion inside the ellipsoid 
D when y = and, for y 7^ 0, to a non-integrable billiard motion inside the perturbed 
eUipsoid Dy. See Figure [TUl 

Applying Theorem|5]to an interior transverse local return map near E"^^, and notic- 
ing that all homoclinic orbits of the billiard flow in Dy are regular orbits, we immedi- 
ately establish: 

Proposition 4. Consider the Hamiltonian flow M-9\ . where y(^;y,e) is a billiard po- 
tential limiting to the billiard flow in Dy (V{q;y,S.) satisfles conditions I-IV for all y 
values). Let the function eo(y) satisfy 

v(8o,y) + m(^*(5(8o,y);eo,y)+M(''(v(£o,y);eo,y)) =o(y) 

and Eo (y) — > ai y ^ 0. Then, for each of the above cases I -4, for sufficiently small 
y > 0, the smooth flow has transverse homoclinic orbits which limit to the billiard's 
transverse homoclinic orbits far allO <E < eo(y)- 
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Figure 10: Perturbation of a billiard flow inside a perturbed ellipsoid family Dy. 

Indeed, for sufficiently small y > equation J48t is valid, and thus the homoclinic 
billiard orbit ^''-hom-y transverse, so the above theorem follows immediately from 
Theorem|5]and the discussion after it. Based on this proposition and Propositions lll3l 
we conclude: 

Proposition 5. Consider the Hamiltonian flow M-9\ . where y(^;Y,8) is a billiard po- 
tential limiting to the billiard flow in Dy (V{q;j,e) satisfles conditions I-IV for all J 
values). Further assume that the potential V{q;j,e) is boundary dominated and is 
given near the boundary of Dy by W{Q;e), so that M-IX holds for the corresponding 
5 values which are specifled bellow. Then, for each of the above cases 1-4, for suffi- 
ciently small Y > 0, the smooth flow has transverse homoclinic orbits which limit to the 
billiard's transverse homoclinic orbits and thus is non-integrable for allO <£ < £0(7). 
where 

• ForW{Q\^) ^e-i : 5 = C»(-elne) and eo{y) =f+'^,x>Q. 

• ForW{Q;z)=e-^ : 5 = 0(V-elne) andtQ{-{) =f+'' ,k > Q. 

• ForW(e;e) = (|)" ; b ^ 0{^+k/^) and =f+^+'^ ,k> Q. 

5 Discussion 

The paper includes three main results: 
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• Theorems 1 1121 deal with the smooth convergence of flows in steep potentials 
to the billiard's flow in the multi-dimensional case. These results, which are 
a natural extension of |35|, provide a powerful theoretical tool for proving the 
persistence of various billiard trajectories in the smooth systems, and vice versa. 
The unavoidable emergence of degenerate tangencies in the higher dimensional 
setting, and the study of corners and regular tangencies (extending L27.LI361 to 
higher dimensions) have yet to be addressed. 

• Theorems 13 141 provide the first order corrections for approximating the smooth 
flows by billiards for regular reflections. Theorem |3] proposes the appropriate 
zeroth order billiard geometry which best approximates the steep billiard and 
a simple formula for computing the first order correction terms, thus allowing 
to study the effect of smoothing. The smooth flow and the billiard flow do not 
match in a boundary layer - the width of it and the time spent in it are specified in 
Theorem^] Propositions 1 113 1 supplv the estimates for the boundary layer width 
and the accuracy of the auxiliary billiard approximation for some typical poten- 
tials (exponential, Gaussian and power-law). All these results are novel for any 
dimension, and propose a new approach for studying problems with relatively 
steep potentials. A plethora of questions regarding the differences between the 
smooth and hard wall systems can now be rigorously analyzed. 

• Theorem 13 and Proposition [S] The above mentioned C' estimates of the error 
terms lead naturally to the persistence Theorem|5l Applying these results to the 
billiards studied in [111, we prove that the motion in steep potentials in vari- 
ous deformed ellipsoids are non-integrable for an open interval of the steepness 
parameter, and we provide a lower bound for this interval length for the above 
mentioned typical potentials. While the analysis of higher dimensional Hamilto- 
nian systems is highly non-trivial, we demonstrate here that some results which 
are obtained for maps may be immediately extended to the smooth steep case. 
We note that the same statement works in the opposite direction. Furthermore, 
one may use the first order corrections developed in Theorem|3land Propositions 
lll3l to study the possible appearance of non-integrabilty due to the introduction 
of smooth potentials. 

As mentioned before, these results may give the impression that the smooth flow 
and the billiard flow are indeed very similar While in this work we emphasize the 
closeness of the two flows, it is important to bear in mind that this is not the case in 
general. This observation applies to the local behavior near solutions which are not 
structurally stable and is especially important when dealing with asymptotic properties 
such as ergodicity, as discussed below. 

Let us first remark about the local behavior First, as in the two-dimensional set- 
tings, we expect that singular orbits or polygons of the billiard give rise to various types 
of orbits in the smooth setting. The larger the dimension of the system, the larger is the 
variety of orbits which may emerge from these singularities. Moreover, in this higher 
dimensional setting, even though our theory implies that regular elliptic or partially- 
elliptic periodic orbits persist, the motion near them (and their stability) may change. 
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Global properties of the phase space are even more sensitive to small changes. If the 
billiard periodic orbit is hyperbolic, while it and its local stable and unstable manifolds 
persist (see for example Theorem|5Jl, their global structure in the smooth case may be 
quite different; First of all, integrability of one of the systems does not imply integrabil- 
ity of the other (for example, it may be possible to use the correction terms computed 
in Section l331 to establish that the smooth flow has separatrix splitting even when the 
billiard is integrable). Second, if the billiard flow has singularities, the global mani- 
folds of a hyperbolic billiard orbit may have discontinuities and singularities whereas 
the global manifolds of the smooth orbit are smooth (see for example 1351 ). 

Finally, the most celebrated global property one is interested in is ergodicity and 
mixing. Indeed, Boltzmann suggested that the gas molecules interacting in a box 
should have a fast decaying correlation function and proposed the analogy of the cor- 
responding dispersing hard balls system. In modern terminology, Boltzmann claimed 
that for sufficiently large systems the hard sphere gases are ergodic and mixing 1211 
and hence so are the real gases. Sinai |29| proved that the dispersion property is suf- 
ficient for proving that the system of two disks on a two-torus is ergodic and mixing, 
and following this fundamental work the study of the dynamics and mixing proper- 
ties in various two-dimensional billiard tables had flourished J5]|6l|38l ?] (the behavior 
of billiards in higher dimensions is much less studied, see 1391 ISl 171 1261 113 1 1 1281 and 
references therein). 

The suspicion that the motion in smooth steep potentials may have a different char- 
acter has been lurking all along. In fact, several works where dedicated to proving that 
in some cases (finite-range axis-symmetric potentials) the motion may be still ergodic 
1291 122I I23I [T51 121 . In 1141 it was shown that when two particles with a finite range 
potential move on a two-dimensional torus stable periodic orbit may emerge. In 1351 
we proved that in the two dimensional case (C smooth potentials, not necessarily fi- 
nite range, not necessarily symmetric), near singular trajectories (tangent trajectories 
or corner trajectories) new islands of stability are born in the smooth flow for arbitrar- 
ily steep potentials. Thus there is a fundamental difference in the ergodic properties 
of hard-wall potentials as compared to smooth potentials. Although these results only 
apply to two-particle systems, they raise the possibility that systems with large num- 
bers of particles interacting by smooth potentials could also be non-ergodic. The tools 
developed here may be useful in studying these possibilities. 
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7 Appendix 



7.1 Proof of TheoremsUandEl 



By Condition I the Hamiltonian flow is C -close to the bilhard flow outside an arbitrar- 
ily small boundary layer. So we will concentrate our attention on the behavior of the 
Hamiltonian flow inside such a layer 

Let the initial conditions correspond to the billiard orbit which hits a boundary 
surface F, at a (non-corner) point qc. By Condition Ila, the surface F,- is given by 
the equation Q{q\Q) = Qi, hence the boundary layer near F,- can be defined as = 
{\Q{q;z) — Qi\ < 5}, where 5 tends to zero sufficiently slowly as e +0. Take e 
sufficiently small. The smooth trajectory enters A^g at some time f,« (5, e) at a point 
q'm(5,e) which is close to the collision point qc with the velocity pi„{5,E) which is 
close to the initial velocity po. See Figure The same trajectory exits from A^g at 
the time foMf(5,e) at a point g'our(5,£) with velocity /?„i,;(5,£). In these settings, the 
theorems are equivalent (r = corresponds to Theorem^] while r > corresponds to 
Theorem^ to proving the following statements: 



lim lim 



^o„r(5,e),f„„r(5,e; 



^,„(5,8),f,„(5,e) 



0, 



(50) 



which guarantees that the trajectory does not travel along the boundary, and (see (|5}) 



lim lim 



Pour (8, e) - p,-„ (5, e) + 2« {qi„) {pi„{5,e),n (qi,, ) ) 



= 0, 



(51) 



where n{q) is the unit inward normal to the level surface of Q at the point q. 

With no loss of generality, assume that Q{q;Q) increases as q leaves D's boundary 
towards D's interior. Choose the coordinates {x,y) so that the hyperplane it is tangent to 
the level surface Q{q;e) — Q{qc',e) and the y-axis is the inward normal to this surface 
at ^ = qc- Hence, the partial derivatives of Q satisfy: 



G.|(,,;e)=0, Qy 



= 1 



(52) 



By (I12> and Condition II, near the boundary the equations of motion have the form: 



dH dH ,, 



(53) 



dH 



-W'iQ;e)Qy. 



(54) 



dpy - - dy 

We start with the C" version of i50\ and i51\ . First, we will prove that given a suffi- 
ciently slowly tending to zero ^(e), if the orbit stays in the boundary layer A^g for all 
t G [tin,tin + ^], then in this time interval 



q{t)^qi„{b,e)+0{^), 



(55) 



P,{t)=p,{ti„{8,e)) + 0{^) 



(56) 



28 




^ + W(6(g(0;e);e) - ^--^^''"^^^'"^^ +wm + Og). (57) 

Note that i55i follows immediately from (I53>-(I54> and the fact that p is uniformly 

bounded by the energy constraint ~ H — W {Q;£,) <H~ j. In fact, qi„ — tends to 

zero as 0(5) for regular trajectories and (9( VS) for non-degenerate tangent trajectories, 
so by assuming that ^(e) is slow enough, we extract from (I55> that 

q{t)=q,: + 0{^). (58) 

Now, from (|53, (|5H}, for t e [f,„(5,e),f,>,(5,e) +^] we have 

e,(<?(0;e) = O©, Qy{q{ty,e) = l+0{^). (59) 

Divide the interval / = [?,„,?,„ +^] into two sets: /< where |W'(2;e)| < 1 and /> where 
|W'(e;£)| > 1. In /< we have = 0{t,) by (|53,(|59}. In />, as |lV'(e;e)| > 1 and 
Qy 7^ 0, we have that py is bounded away from zero, so in i53\ we can divide p^ by py: 

dpx _ & 

dpy Qy 

It follows that the change in p^ on / can be estimated from above as (the contri- 

bution from /<) plus 0{t,) times the total variation in py. Thus, in order to prove i56\ . 
it is enough to show that the the total variation in py on / is uniformly bounded. Recall 
that Py is uniformly bounded (\py\ < 1 from the energy constraint) and monotone (as 
^'(2) < and Qy > 0, we have py > 0, see i54i ) everywhere on /, so its total variation 
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is uniformly bounded indeed. Thus, i56\ is proven. The approximate conservation law 

i5T\ follows now from (|56j and the conservation of H = ^ + ^ + W{Q{q;E);E). 

Finally, we prove that Tg, the time the trajectory spends in the boundary layer A^g, 
tends to zero as £ ^ 0. This step completes the proof of Theorem ^ by plugging 
the time Tg ^ instead of t, in the right-hand sides of ( I55t .( l56> .( l57t . we immediately 
obtain the C''- version of iSOi and i5H . 

Let us start with the non-tangent case, i.e. with the trajectories such that Py{tm) is 
bounded away from zero. From Condition III it follows that the value of W,„ = Wuut = 
W{Q = 5;£) vanishes as e — > +0. Hence, by ( I57> the momentum Py{t) stays bounded 
away from zero as long as the potential W{Q^,^) remains small. Choose some small v, 
and divide A^g into two parts A^< := {W : W{Q;e) < v} and N> = {W : W{Q;e) > v}. 
First, the trajectory enters N^. Since the value of -jjQiq) — PxQx + PyQy is negative 
and bounded away from zero in (because 2a is small, and py and Qy are non-zero), 
the trajectory must reach the inner part by a time proportional to the width of A^<, 
which is (9(5). Also, we can conclude that if the trajectory leaves A^> after some time 
f>, it must have py > and, arguing as above, we obtain that tout — tin — 0{5) +f>. Let 
us show that f> ^ as e ^ +0. Using i54i . the fact that the total variation of py is 
bounded, and Condition IV, we obtain 

C 

k>l<^ , =Cmax\(l'(W-e)\^0 ase->+0. 

' - mmwJW'(e;e)| iv> ^ ^' 

So, in the non-tangent case, the collision time is (9(5 + f>), i.e. it tends to zero indeed. 

This result holds true for py ,„ bounded away from zero, and it remains vahd for 
Pyjn tending to zero sufficiently slowly. Hence, we are left with the case where p,, ,„ 
tends to zero as e (the case of nearly tangent trajectories). Inside A^g, since W is 
monotone by ([l9j, we have W(g;e) > Wi„ = W(5;e). Therefore, by (EJ, Py{t) stays 
small unless the trajectory leaves A'g or t — tj„ becomes larger than a certain bounded 
away from zero value. From \56\ it follows then that px{t) remains bounded away from 
zero. By (I53t.(l54l. 

Q-=^Q{q(t)\^)^QxPx + Q.Py 

at 

so Q is small, yet 



— Q{q(t)-t) = plQxxPx+2Q„PxPy + Qyyp]. " W {Q-€){QI + Q]) . 

For a non-degenerate tangency, p^QxxPx is positive and bounded away from zero. 
Therefore, as py is small and W'{Q;e) is negative, we obtain that ^Q{q{t)\^) is posi- 
tive and bounded away from zero for a bounded away from zero interval of time (start- 
ing with tin). It follows that 

Q{q{t);^)>Q{q,„;^) + Q{tin){t-t,n) + C{t-t,„f (60) 

on this interval, for some constant C > 0. We see from ( I60> . that the trajectory has to 
leave the boundary layer A^g = {\Q{q;£) — Qi\ < 5 = \Q{qi„;t) — 2,- 1} in a time of order 
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0{Q{tin)) = 0{Q^{qi„))+0{py^in) = 0{qi„-qc)+0{pyjn)- As qi„-qc = 0{V5) for 
a non-degenerate tangency, we see that the time the nearly-tangent orbit may spend 
in the boundary layer is 0{V5 + Pyj,,), i-C- in this case it tends to zero as well. This 
completes the proof of Theorem^ 

Now we prove Theorem|2l- the C -convergence for the non-tangent case. Again, 
divide into A^< and A^> for a small v and consider the limit lims^o limy^o limg^+o- 
As we have shown above, 2 7^ in A^<, thus we can divide the equations of motion 
O, (|53 by Q: 



dq p 

dQ " QxPx+PyQy 

dp I, ^ VQ 

— = -W'(e;e) , 

dQ QxPx+PyQy 



(61) 



dt 



1 



dQ QxPx+PyQy 

Equations ( 16 U can be rewritten in an integral form: 

rQi 



q{Q2)-q[Qx)^ / F,iq,p)dQ, 



piQ2)-piQi) 



W{Q2 



Fp{q,p)dWiQ), 



(62) 



t{Q2)-t{Qi)^ Ftiq,p)dQ, 

where Fy , Fp and Ft denote some functions of {q,p) which are uniformly bounded along 
with all derivatives. In A^<, the change in Q is bounded by 5 and the change in W is 
bounded by v. Hence, the integrals on the right-hand side are small. Applying the 
successive approximation method, we obtain that the Poincare map (the solution to 
i62\ ) from Q = Qi to Q = Q2 limits to the identity map (along with all derivatives with 
respect to initial conditions) as 5,V ^ 0. It follows that in order to prove the theorem, 
i.e. to prove (Bol.dsTl. we need to prove 



Um lim 



0, 



(63) 



and 



lim lim 



Pout Pin 

+ 2n{qi„){pi„,n{qi„)) 



(64) 



where {qm, Pin,tm) and {q our, Pout, tout) correspond now to the intersections of the or- 
bit with the cross-section W{Q{q,e),e) = v. By Condition IV, as e — > the function 
Ci{W;e) tends to zero uniformly along with all its derivatives in the region v <W < H 
for any v bounded away from zero. Therefore, the same holds true for a sufficiently 
slowly tending to zero V and W'{Q;e) — (Q,'(VK;e))^' is bounded away from zero in 
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the region A'> . Hence, by J54> . the derivative py is bounded away from zero as well. 



Therefore, we can divide the equations of motion (I53>.(l53ll by 

;^ = -cz'(W;s)-f, ^—ci'iW;e)^, ^ = |, (65) 

dpy Qy dpy Qy Clpy Qy 



W^H-^p\ (66) 



where 



Condition IV implies that the C -limit as £ ^ of i65i is 

^^0, (67) 

dpy dpy Qy 

Since the change in py is finite and the functions on the right-hand side of \65\ are all 
bounded, the solution of this system is the C^-limit of the solution of \b5l . From 
we obtain that in the limit e ^ {qin,ti„) — {qout,tout), so ( I63> is proved. Second, we 
obtain from (|^} that 

{Px.oiit ~ Px.in)Qy{qm\^) = [Py.oiit - Py,in)Qx{qin\^) 

in the limit e ^ 0, which, in the coordinate independent vector notation (see e.g. I34> . 
and by using (<7,>,,f,>,) = {q out, tout), amounts to the correct reflection law. 

7.2 Picard iteration for equations with small right-hand side. 

Before we proceed to the proof of Lemmas ^ and |3] we recall the main tool of their 
proofs - the Picard iteration scheme for equations with small right hand side. Consider 
the differential equation 

v = i|/(v,^,f,e) (68) 

where \(/ is a C -smooth function of v and fi, continuous with respect to t and e. Assume 
that for t e [0,L(£)] and bounded {v,/j) we have a function J{e) such that 7(e)L(e) 
and 

||v||c'-</(e). (69) 
Then, according to the contraction mapping principle, the Picard iterations v„ where 



v„+i{t) ^vo+ [ \\f{v„{s),ij,s,e)ds 
Jo 



(70) 



converge to the solution of ( I68> starting at f = with initial condition v(0) = vq on the 
interval f G [0,L(£)], in the C^-norm as a function of vq and /j: 

v„{t;vo,ij) —f , v{t\VQ,n) = V0+ / ^{v{s\VQ,n),n,s,€)ds. 

Jo 

One can show by induction that ||y„(f) — vo||^r = (9(L(e)7(e)) uniformly for all n. Then 
it follows that 

v{t-vo,ij) = vo + 0^,iL{e)J{e)). (71) 
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It is easy to show that 

v(f;vo,A<) =v„(f;yo,A') + 0^.-,((i(e)-/(e)r+') (72) 

(such kind of estimates are, in fact, a standard tool in the averaging theory). In order 
to prove i72\ . we will use induction in n. At n = we have even better result than j72> 
(see (17 H). Now note that 

v{t) ~ v„+i{t) = [ {\\i{v{s),ij,s,e)~\\t{v„{s),ij,s,e)ds 



""i) \Io 'VvMs)+z{v{s)-v„{s)),p,s,e)dzj ■iv{s)-v„{s))ds. 
It follows immediately that 

v„+i||^,._, = 0(L(e)|K||^,._, ) • 0(|lv- v„||^.,._, ) = 0(L(8)7(8)) • ||v- y„||^.,,_, , 
and ( I72> indeed holds true by induction. 

7.3 Proof of Lemma [l] 

The "free flight" (the motion inside D^) is composed of motion in D?,^ (the region 
outside of N^) and the motion in the layer A^< ~ D^\D^^^j. We show that in each of these 
regions the equations may be brought to the form ( I68> ,(I69>. We will first consider the 
flight inside Recall that the equations of motion for the smooth orbit are 

(73) 

Let us make the following change of coordinates 

q{t):^q{t)~p{t)t (74) 

Then (I73> takes the form 

p = -VV{q + pt;e) 

with initial data {q{Q),p{Q)) — {qo,po)- Since the time spent in D^^^ must be finite as it 
is C — close to the billiard's travel time in Df^^ which is finite here, and using ( I28> . we 
have 

Thus, system J75> does satisfy i69\ with L = (9(1), J = Olm^''^). It follows then from 
(ED that 

p(f)=po + 0^. (76) 

Furthermore, by applying « = 1 Picard iteration (I70> . we obtain from M2\ the following 
estimate for p{t): 



p{t) =P0- ^^(?o + pos\e)ds + 0^.,._, ((mW)2). 



(77) 
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By integrating the equation p,we also obtain from J76t that 

q{t)=^qo + pQt + O^,- (m*'') ) . (78) 

Next, we show that the equations in the layer A^< = {W : W{Q;z) < v} can be brought 
to the form ( I68t .(l69ll as well. Recall (see the proof of Theorem|2ji that Q = {VQ,p) is 
bounded away from zero in A^<, hence 2 can be taken as a new independent variable (it 
changes in the interval V[ < Q — Qi < 5). Now the time t is considered as a function of 
Q and of the initial conditions {q{ti,),p{ti,) (where f§ is the moment the trajectory enters 
A^<). Recall that we showed in the proof of Theorem |2] that f is a smooth function of 
the initial conditions, with all the derivatives bounded. So, in A^<, we rewrite i75\ as 

dQ - ^ y^^^) {VQ{ci+pt:e).py 
dp_ _ MrHr,.,-\ V6(^+p?;e) 
dQ - lt^'tj(ve(^+pr;e),p)- 

As W is a monotone function of Q (i.e. W'(2;e) ^ 0), we can take W as a new inde- 
pendent variable, so the equations of motion will take the form 

dq_ ^ NQ[q+pV,Z) 

dW {VQ{q+/n-£).p)' frjQ. 

dp _ VQ{q+pt;E) ^'^> 

dW ~ {VQ(q+pt;e),p) ' 

Since all the derivatives of f with respect to the initial conditions are bounded, we may 
consider il9i as the system of type ( I68> ,(l69ll with J = (9(1), and L — 0(v) (recall that 
the value of W changes monotonically from Wo = W (5, £) to v). Thus, by applying one 
Picard iteration ilOt . we obtain from j72t that 

' ' iw(8,£) (Ve(^(Wo)+/?(Wo)f;e),/?(Wo)) ^--'^ ' 

From M\\ we also obtain 

p{W)=p{Wo)+0^,(v). 

Note that O^^ ^ (v^) and (9^, (v) refer here to the derivatives (with respect to the initial 
conditions) of p at constant W or, equivalently, at constant Q. Returning to the original 
time variable, these equations yield 



pit) = p{t^) + O^, (v) = p{t^) - f Vy (^(fg) +;,(f5) {s - f8);e)di + Ocr-i (v^). 

•''6 

Using expressions ( I76> . (I77> for /^(fs) and ( I78> and qit^), we finally obtain 

p{t) = po + O^., (V + m(')) ^PO- ^Viqo + pos;e)ds + Ocr-i {{m^''^ +vf) (80) 

for all t such that q{t) E D^, in complete agreement with the claim of the lemma (as we 
mentioned, the O , j (•) and (9^, (•) terms refer to the derivatives at constant Q). The 
corresponding expression for q{t) (see i30i ) is obtained by integrating the equation 
q = P- The expression (13 1> for the flight time T is immediately found from the relation 
W{Q{q{z);E);E) = V or, equivalently, Q{q{t);e) = g, +r| (recall that Q is bounded 
away from zero in the layer A^<). 
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7.4 Proof of Lemma 121 

Here we compute the reflection map 7?^ : {qi„,pm) {qout,Pout) defined by the smooth 
trajectories within the most inner layer A^> : {W > v}. We put the origin of the coordi- 
nate system at the point qi„ (corresponding to q at Figure|5} and rotate the axes with e 
so that the y-axis will coincide with the inward normal to the surface Q{q; e) = Q{qin\ e) 
at the point qi„ (corresponds to n{q) at Figure|5}, the x-coordinates will correspond to 
the tangent directions. It is easy to see that in the notations of Lemma|3]we have (the 
explicit dependence on e is suppressed for brevity) 

K{qin) = Qxx{qin)/Qy{qm), Kqin) = Qxy{qin) / Qy{qm) ■ (81) 

As we have shown in the proof of Theorem |2] ^ is bounded away from zero in A^>; 
Hence, we may use Py as the new independent variable (see j65t). In order to bring the 
equations of motion to the required form with the small right hand side, we make the 
additional transformation 

Px ^P^Px- TTTTPy ■ (82) 
Note that Qx{qi„\s) ~ 0, hence (see dsTl ) 

P = Px-K{qi„){x-Xin)py - l{qm){y-ym)Py + 0{{q - qinf). (83) 
In particular 

p{tin)=PxM- (84) 

After the transformation, equations (165 > take the form 



;^ = -(1'(t-tp')7^, (86) 

dpy 2 2 Qy 

dp ^ 2\d ( Qx\ P 



dpy 2 2 dq \ Qy J Qy 

Since Q,'(VK;e) is small in the inner layer, these equations belong to the class (I68ll.(l69>. 
with J ^ C>(m('")) (see ^) and L = 0(1) (the change in Py is bounded by the energy 
constraint). Thus, by Mil , we obtain (see i&AV ) 

{q,t,p) = {qi„ , ti„ , p^^i„ ) + O^, (m''') ) . (88) 

Recall that 'W{qout) — W{qin)- Therefore, by energy conservation, 

Px.in ~l~ Py.in Px.out ~^ Py,out i (89) 

SO ( I88t implies 

Py,aut = -Pyjn + O^, [M^'^ ) . (90) 
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BvdSSt. and by using Qx{qm,^) =0, equations dSSt may be written up to ((m'''')^)- 
terms as 

apy 2 Qy\qin) 

(92) 

^ = a'(5(i-pi,.-p;))(ifte,,)p.j«+'(9/Jp,))g^- (93) 

Now, by applying to equations \%5\ the estimate (I72> with n = 1 (one Picard iteration), 
we can restore from (19 1> all the formulas of lemma |3] (we use (I83> to restore px from 
p, and use ( I89> to determine Py.out', note also that, up to (9(M('^')-terms, the interval of 
integration is symmetric by virtue of ( I88> . so the integrals of odd functions of py in the 
right-hand-sides of (|^ are C>((m(''))2)). 
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